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We calculate the dc conductivity cr in a pseudogapped high Tc s perconductor within a precursor 
superconductivity theory which is consistent with gauge invariance. Our results contain physical 
effects beyond those identified previously. Rather than presuming that lifetime effects dominate the 
T dependence of transport, here we show (consistent with growing experimental support) that the 
temperature dependence of the effective carrier number is a natural consequence of the pseudogap, 
and demonstrate reasonable agreement with dc transport in the underdoped cuprates. 
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Understanding the T dependence of the dc resistiv- 
ity, particularly near optimal doping, was one of the 
first puzzles posed by the high temperature supercon- 
ductors. A number of different models based on spin- 
charge separation to marginal Fermi liquid phenomenol- 
ogy were invoked to explain the unusual power laws ob- 
served. With our increased understanding of the pseu- 
dogap phase (which extends over most of the phase di- 
agram, including optimal doping), it has become clear 
that earlier theories must be modified to accomodate the 
pseudogap [1] and related Fermi arc effects. Moreover, 
the nature of this pseudogap is currently under heated 
debate as to whether it derives from precursor and/or 
fluctuation superconductivity [2] or an alternative order 
parameter [3]. It is therefore essential to build a sys- 
tematic and properly conserving transport theory within 
each of these schools and their variants. 

In this paper we use a pre-formed pair, sum rule- com- 
patible approach to address the dc conductivity cr(r) in 
the cuprates. Our goals are to determine (i) the ef- 
fects of a pseudogap on a(T), (ii) the experimentally 
deduced temperature dependence of the lifetime, (iii) 
how to understand superconducting coherence effects in 
a{u! 0,T) [4] from above to below Tc, (iv) the role of 
the Fermi arcs as well as (v) what is the source of the 
"bad metallicity" [5] in the cuprates. Our approach is 
based on the notion that the attraction is stronger than 
in BCS theory; this leads to finite center of mass pair 
excitations of the normal and the condensed state. The 
formal machinery used here for addressing transport is 
well established [6-8]. These pair fluctuations are a nat- 
ural consequence of the short coherence length, as well 
as generally high Tc. They are to be distinguished from 
phase fluctuations; as a mean field approach we do not 
include superconducting fluctuations which dominate in 



the critical region. 

For illustrative purposes, it is helpful, first, to present 
our expression for the dc conductivity for s-wave jellium, 
and thereby illustrate what is the role of pre-formed or 
non-condensed pairs. We find 

Here E'k is the usual BCS dispersion, E'k = v^k+^k 
where the excitation gap consists of two contributions 
from non-condensed (pg) and condensed (sc) pairs: via 

= A^^ ,^ -I- A^^^, f{Ep) is the Fermi function and 
is the effective lifetime. The gap Ak remains relatively 
T-independent, even below Tc, as observed, because of 
the conversion of non-condensed (Apg.k) to condensed 
(Asck) pairs as the temperature is lowered. This trans- 
port equation is for the weak scattering limit and when 
Apg = it reduces to strict BCS theory, where below Tc, 
we consider cj — > 0+. 

Importantly, this scheme is associated with the stan- 
dard [9-11] pseudogap self energy, which we derived even 
earlier within our microscopic formalism [12] 

Ep,(if) = -z7+-^f^ (2) 

Here 7 represents a damping related to the inter- 
conversion of pairs and fermions. (Rather than intro- 
ducing two different 7's in the first and second term, 
we minimize the number of parameters and take them 
equal.) 

Physically, Eq. (1) implies that the dc conductivity is 
affected by the pseudogap or pre-formed pairs, via a re- 
duction as well as a temperature dependence of effective 
carrier number (n/m(T))cff. The first of these arises 
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from the Fermi function derivative, which reduces the 
number of carriers through a gap effect. The second of 
these appears via the pre-factor — A^^ which reflects the 
fact that when fcrmions arc tied up into pairs the car- 
rier number is decreased. Indeed, there is an increasing 
experimental awareness [13, 14], that "the do resistivity 
of cupratcs pdc{T) is governed not only by the relaxation 
processes but also by temperature-dependent numbers of 
carriers" . 

Formalism A central question here is how to incorpo- 
rate the widely used Eq.(2) into a consistent treatment 
of transport. Based on Eq.(2), one can write for the full 
Green's function including condensed pair (sc) effects 
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Note that there can be no finite lifetime effects associated 
with the condensed pairs. 

Within a BCS-like formulation, transport properties 
will involve terms of the form Fsc.k which represent the 
usual Gor'kov functions as a product of one dressed (G) 
and one bare (Go) Green's function (GGq), 
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To address transport we consider the EM kernel J = 
-Va, with K{Q) = e2(n/m).. +P{Q), where the 

paramagnetic contribution, given by P{Q), is associated 
with the normal current from fermionic and bosonic ex- 
citations. 

The diamagnctic current (with electronic dispersion 
Ck), is iU/m)dia = 2^^(52^k/5k5k)Gx. We integrate 
this expression by parts and use the generalized Ward 
identity to obtain [7] an alternate form (above T^) 
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In applying the Ward identity to arrive at Eq. (5), we 
are considering a t-matrix approach in which the self en- 
ergy can be written as T.pg{K) = '^Qtpg{Q)Go{Q — K), 
with tpgiQ) = ([/-I -I- Xpg{Q)y^ where XpgiQ) = 
Y^k^kGq^q-k- Here U is the strength of an attractive 
interaction which is unspecified. The presence of GGq 
in the non-condensed pair t-matrix, tpg follows from the 
GGo form of the Gor'kov function (Eq. (4)). We de- 
fine K = (k, {Q = (q, ifim)) where iw„ {i^m) is a 
fermionic (bosonic) Matsubara frequency. 

Eq.(5) is important because it has cast the diamag- 
netic response in the form of a two particle correlation 
function. That there is no Meissner effect in the nor- 
mal state is related to a precise cancellation between the 
diamagnetic and paramagnetic terms. Noting P(0) = 
— e^(n/m)dia) we may extend -P(O) to finite Q and infer 



that in the normal state 

t^Ck-l-q/2 5^k-|-q/2 

K 
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p 

One can alternatively [7, 8] introduce the Aslamazov- 
Larkin and Maki-Thompson diagrams to arrive at the 
above equation [9]. The condensate contribution is sim- 
ilarly associated with a t-matrix. T,{K) = T,sc{K) + 
^pg{K) and 11{K) = Y.Qt{Q)Go{Q - K). where t{Q) = 
tsrXQ) + tpgiQ). Since tsc{Q) = -S{Q)Al^ jT, this 
yields the superconducting contribution in Eq. (3) 
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Then, in the same spirit as our previous [12] deriva- 
tion of Eq.(2) we take note of the fact that up to 
T w r*/2, where T* is the pairing onset temperature 
(ie, where the pair chemical potential is small), tpg(Q) 
is strongly peaked at small P. This leads us to identify 
Apg — — J2 tpgiQ) and clarifies the parameters of Eq. (2). 
We rewrite Eq. (6), now including the usual condensate 
contribution as 



PiQ) ^ 2e'J2 
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+ Fsc,kFsc,K+Q — Fpg^xFpg^K+Q 



with Fpg,K = - . ^Pf'', . Gk 



(8) 



(9) 



From Eq.(8) and Re(Tp°''"(n) = (^^^L^^^), the para- 
magnetic contribution to the dc conductivity is 



-Fpg^KFpg^Ky:i + Fsc,KFsc,K.\Q^ 



in^->-o+ 



.(10) 



In previous work in the literature [11, 15] only the first 
term involving GG was included, which was recognized 
as inadequate [11]. Note also that Eq (1) follows directly 
from Eq. (10) in the limit of small 7 « t~^. For consis- 
tency we rewrite Eq. (5), also adding in the usual BCS 
condensate terms 
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We are now in a position to demonstrate compat- 
ibility with the important conductivity sum rule 
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Figure 1: (Color online) Conductivity in an overdoped sample 
(a) The effective carrier number. The inset shows the effective 
carrier number for 7/A(rc) = 0. (b) a{T/Tc). (c) Different 
contributions to 'Rea{T/Tc): the contributions due to Asc 
and Apg are labeled and the total value of a{T/Tc) is the 
unlabeled region. Panels (b) and (d) compare favorably with 
experiment [4, 16]. 



dnRea{^) = {l/2)e^{nxx/'m)dia- Note that 
Rea{n) = -ImPxx{n)/n + Tr5{n)[RePxx{n) + 
(^icaj/'Ti)^^^]. Integrating the first term over fre- 
quency we find = —TrPxx{0), while the second (delta 
function) term yields a term +ttPxx{0), which leaves 
only the diamagnetic contribution and yields the desired 
sum rule. Note that this sum rule is intimately connected 
to the absence (above T^) and the presence (below Tc) 
of a Meissner effect. Importantly, since {~^)d^^ can be 
viewed as essentially independent of temperature, when 
there are approximations in evaluating the transport 
diagrams, it is appropriate to evaluate the chemical 
potential /i based on the T- independence in Eq(ll). 

Numerical Results For our numerical calculations 
we have essentially two fitted parameters T* and Tc, rep- 
resenting hole doping, and these in turn constrain the 
attraction U. We obtain n and Ak(T) from the mean 
field equations where the mean field transition temper- 
ature is T*. Rather than using a more numerically in- 
tensive approach [9] we simplify the decomposition of pg 
and sc contributions so that \g,kiT) = {T/T^f/'^A]^{T) 
for T < and for T > T^, ' Apg^^iT) = At- This 
parametrization of the gaps is very close to the full 
numerical solution [9]. The values for 7 were chosen 
to give rough agreement with the observed Fermi arc 
lengths. Although it is not particularly critical, we take 
the form 7 oc {T/Tcf for T < T^. This is based on 
fitting the temperature dependence of the spectral func- 
tion peak width to the measured quasi-particle scatter- 
ing rate [17]. We define the effective carrier number 
{n/m{T))oB = j{T)a(T). It should not be surprising 
that the effective carrier number increases with temper- 



Figure 2: (Color online)Transport in an underdoped sample, 
(a) The effective carrier number. The inset shows the effective 
carrier number for "f/A{Tc) — 0.(b) g{T/Tc). (c) Various 
contributions to 'RecriT /Tc). Corresponding (d) resistivity. 
Panels (b) and (d) compare favorably with experiment [4, 16]. 



4 I — I — I — I — I — I — I — I — I — I — I — I — I — r 




Figure 3: (Color online) Doping dependence of resistivity (be- 
low T*) and (inset) ARPES gaps showing Fermi arcs. The 
normalization po is defined as p{Tc) for the highest doping 
case (7/A(rc) = 0.67). The change in doping is associated 
with different ratios of j/AiTc) = 0.07, 0.20, 0.35, 0.67 for in- 
creased doping. 

ature because of the presence of an excitation gap and 
related non-condensed pairs. This T-dependent carrier 
number is precisely what is found in a strict BCS super- 
conductor below Tc- This contribution leads to a non- 
metallic tendency with a increasing with T above Tc. 

Experimentally [1], one finds an increasing resistivity 
with T and a quasi-linear T dependence. Convention- 
ally this is explained by assuming that (ri/m(T))og is 
T- independent, but that the inverse lifetime is linear in 
temperature. We have just shown that due to the pseu- 
dogap, the carrier number necessarily increases with T. 
This suggests that 7(T) must be a higher power than 
linear. It is not implausible to assume that the gapless 
Fermi arcs lead to a more conventional, Fermi liquid be- 
havior; thus we choose 7 oc (T/Tc)^. 

In Fig. 1 we plot results for a prototypical overdoped 
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system with T* /T^ = 1.4 and 7/A(Tc) = 0.67. In the up- 
per left inset we also plot {n/m{T))es in the limit that 
7 — >^ 0, so that there is no distinction between condensed 
and non-condcnscd pairs and Fermi arcs arc not present. 
The lowest T considered is T/Tc = 0.4 where the numer- 
ics are well controlled and impurity effects are less impor- 
tant. The highest temperatures considered correspond to 
T*, above which the simple form for the fermionic self en- 
ergy of Eq. (2) is no longer valid. It should be noted that 
both cr and (n/?Ti(T))off vanish at strictly zero tempera- 
ture, as we do not include a self consistent treatment of 
d-wave localization [18]. Plotted in Fig. 1(a) is the cor- 
responding effective carrier number (n/m(T))off. From 
this figure one can see that the carrier number begins to 
rise once T exceeds Tc- Because the difference between 
the inset and main figure are not dramatic, we see that 
the arcs are not crucial for this aspect of transport. In 
Fig. 1(b) we plot Rea-(T/Tc), which exhibits a maximiun 
below Tc, as has been been observed experimentally [4]. 
This is associated with the onset of order via Age, which 
increases a (Eq. (10)). Above Tc, we see that the con- 
ductivity of the normal state is appropriately metallic, 
but suppressed by the excitation gap. In Fig. 1(c) we 
plot and label the Age and Apg components of Recr aris- 
ing from the three different contributions on the right 
in Eq. (10) as well as the total dc conductivity. The 
(pg) contribution from the non-condensed pairs lowers 
the conductivity because the presence of non-condcnscd 
pairs means fewer fermions are available for dc transport. 
The resistivity p, shown in Fig. 1(d) has a nearly linear 
temperature dependence, but rises faster than T near T* 
because the effective carrier number is nearing its normal 
state value. 

The counterpart results for an undcrdopcd sample with 
a much larger gap and 7/A(Tc) = 0.07 and T*/Tc = 4, 
are shown in Fig. 2. The effective carrier number is plot- 
ted in Fig. 2(a), and the inset shows the results for 7 = 0. 
Above Tc, the curvature of (n/m(T))eff changes to con- 
vex bowing, whereas it is concave in the overdoped case. 
This is to be expected as this larger value for A{Tc) sup- 
presses the carrier number. The conductivity plotted in 
Fig.2(b) is similar to its experimental counterpart and 
the difference between this and the overdoped case is 
made apparent in Fig. 2(c). One sees that the contri- 
bution from non-condensed pairs is much larger at lower 
doping due to the larger pseudogap. The larger depres- 
sion in (n/TO(T))eff in turn leads to a concave bowing in 
the resistivity, seen in Fig.2(d). 

Figure 3 presents a plot of normal state resistivities 
and (inset) spectral gaps associated with the fermionic 
spectral function for dopings that interpolate between 
the undcrdopcd and overdoped cases. The resistivities 
are normalized by the value po = p(Tc) for the case 
of highest doping. We may characterize each curve, 
in order of increasing doping, by the ratio 7/A(Tc) = 
0.07,0.20,0.35,0.67, which lead to progressively larger 



arcs. The size of the resistivities decreases as one in- 
creases the hole concentration and this reflects the change 
in gap size and effective carrier number. One can see 
from the figure that there is a change in the nearly lin- 
ear slope with increased doping from concave to convex 
bowing, which may have been seen experimentally [16]. 
The spectral gaps at Tc, displayed in the inset of Fig. 3 
show the general experimental trend [1, 11] where the arc 
length increases with doping. 

Conclusions Our main observation of T-dependence 
in the carrier number (and its implications for dc trans- 
port) seems quite natural and general, although the 
community focus has been on the T dependence of the 
scattering time. The present pre-formed pair approach 
is importantly demonstrably consistent with conserva- 
tion laws and extendable below Tc (merging into a BCS 
ground state). By contrast, other papers in the literature 
[2, 11, 15] are restricted to T > Tc and, unlike Ref. 3 do 
not establish sum rule consistency. The collapse of the 
Fermi arcs and the so-called two-gap physics (through 
Eq. (3)) are also evident [19] here, as is the anomalously 
low conductivity, or "bad- metal" behavior [5] . The (gap- 
less) Fermi arcs appear to be more important for the 
relaxation time than for the carrier number. 
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